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Abstract 

In the present work we show, analytically and numerically, that the variance of many-particle 
operators and their uncertainty product for an out-of-equilibrium Bose-Einstein condensate (BEC) 
can deviate from the outcome of the time-dependent Gross-Pitaevskii dynamics, even in the limit 
of infinite number of particles and at constant interaction parameter when the system becomes 
100% condensed. We demonstrate our finding on the dynamics of the center-of-mass position- 
momentum uncertainty product of a freely expanding as well as of a trapped BEC. This time- 
dependent many-body phenomenon is explained by the existence of time-dependent correlations 
which manifest themselves in the system’s reduced two-body density matrix used to evaluate the 
uncertainty product. Our work demonstrates that one has to use a many-body propagation theory 
to describe an out-of-equilibrium BEC, even in the inhnite particle limit. 
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I. INTRODUCTION 


The out-of-equilibrium dynamics of a quantum system is described by the time-dependent 
Schrodinger equation. All physical information on the evolving quantum system can thus 
be obtained from its time-dependent wavefunction by applying various operators and calcu¬ 
lating expectation values. The variance of an operator quantihes to what extent the system 
nnder investigation is in an eigenstate or a snperposition of eigenstates of the operator. In 
this sense it dictates the quantum resolution by which the operator could be measured. The 
product of the variances of two operators dehnes an uncertainty product. The uncertainty 
prodnct qnantihes to what extent two operators can be mntually measnred. As such, it is a 
fundamental concept in qnantnm mechanics. A famous example is the position-momentum 
uncertainty product of a single quantum particle which is analyzed in qnantnm mechanics 
textbooks for both the static and dynamic cases, see, e.g., [1]. 

Over the past two decades, since they were first experimentally realized [2-6], Bose- 
Einstein condensates (BECs) made of ultracold trapped bosonic atoms have become a pop¬ 
ular ground to stndy interacting qnantnm systems, see the reviews [7-9] and books [10-12]. 
There has been an intense theoretical interest in BECs, and ample studies have been made 
to describe their static and particularly dynamic properties using Gross-Pitaevskii, mean- 
held theory. The time-dependent Gross-Pitaevskii equation governs a mean-held theory 
which assumes that each and every boson is described by one and the same time-dependent 
one-particle fnnction throughout the evolution of the BEG in time. 

The general paradigm is that Gross-Pitaevskii theory properly describes the gronnd state 
as well as the out-of-equilibrium dynamics of BECs in the limit of large particle nnmber. 
To be explicit, Lieb, Seiringer, and Yngvason have rigorously proven for trapped BECs 
with two-body repulsive interaction, in the limit of inhnite particle nnmber and at constant 
interaction parameter, that the ground-state energy and density of the condensate converge 
to those obtained by minimizing the Gross-Pitaevskii energy functional [13]. Thereafter, 
Lieb and Seiringer proved in the same limit that the ground state is 100% condensed [14]. 
In the case of out-of-eqnilibrium dynamics, Erdos, Schlein, and Yau have rigorously proven 
that an expanding initially-trapped BEG, after the trap is released, still exhibits 100% 
condensation [15]. Furthermore, the condensate density evolves according to that predicted 
by the time-dependent Gross-Pitaevskii equation. 
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In a previous work [16] we analyzed the ground state of a trapped BEC and demonstrated 
that, even in the inhnite particle limit when the BEC is 100% condensed, the variance of a 
many-particle operator can substantially differ from that predicted by the Gross-Pitaevskii 
theory. The existence of many-body effects beyond those predicted by the mean-held Gross- 
Pitaevskii theory stems from the necessity of performing the inhnite particle limit only 
after the quantum mechanical observalbe is evaluated and not prior to its evaluation. This 
is essential since otherwise any trace of many-body correlations is washed-out before the 
quantum mechanical observable can be evaluated. This has been explained in length both 
analytically and numerically in Ref. [16]. 

In the present work we generalize our previous result to the dynamics of an out-of- 
equilibrium BEC. Dynamics is generally more intricate than statics, and involves (many) 
excitations. We show, analytically and numerically, that the evolution in time of the un¬ 
certainty product of two operators can deviate from that of the Gross-Pitaevskii dynamics, 
even in the inhnite particle limit. We explicitly demonstrate this deviation for the center- 
of-mass position-momentum uncertainty product of a freely expanding BEC as well as to 
the dynamics of a trapped BEC. Our work advocates that one has to use a many-body 
propagation theory to describe the out-of-equilibrium dynamics of BECs, even in the limit 
of inhnite number of particles when the system becomes 100% condensed. 

The structure of the paper is as follows. In Sec. II we analyze the uncertainty product 
of a freely expanding, initially-trapped BEC. In Sec. Ill we present a general theory for the 
many-body uncertainty product in trapped BECs, and show its importance for the case of 
the bosonic-Josephson-junction system. Concluding remarks are put forward in Sec. IV. 


II. THE UNCERTAINTY PRODUCT OF A FREELY EXPANDING BOSE-EINSTEIN 
CONDENSATE 

Even for the simplest case of an expansion from an harmonic trap, one expects in the limit 
of inhnite number of particles and at constant interaction parameter the Gross-Pitaevskii 
and many-body dynamics to coincide. Whereas following Ref. [15] this indeed holds for 
the density of the BEC, as we shall see below this no longer is the case for the uncertainty 
product. The expansion of a BEC from an harmonic trap has amply been discussed in 
the literature at the mean-held level, see., e.g.. Refs. [18-20]. Furthermore, the uncertainty 
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product of an expanding BEC can be analyzed rather straightforwardly at the many-body 
level due to the separability of the center-of-mass coordinate. Yet, as far as we know, 
contrasting at the inhnite particle limit the exact time-dependent uncertainty product and 
the Gross-Pitaevskii one has not been made, which is the purpose of this section. For the 
sake of completeness, we give a short derivation of the relevant many-body quantities below. 

Consider the many-body Hamiltonian of N interacting bosons in a three-dimensional trap 
Vrir), 


N r 




i=i 


1 

' 2 ^ 


+ VriTj) 


^ AolY(rj - Ffc) 


j<k 


N 


+ Y^T{rj). 
i=i 


( 1 ) 


Here, h = m = 1, and W{yi — r 2 ) is the boson-boson interaction with Aq its strength. Here 
we only consider the case of repulsive interaction, i.e., Aq > 0. The Cartesian components 
are Vj = {xj,yj,Zj), ^ ■ The Hamiltonian can in principle be time 

dependent, although we will not exploit this option explicitly below. 

The system is typically, but not necessarily, prepared in the ground state of the trap 
Vrir), 

Hriri, ... ,rAr)<h(ri,..., ftv) = E<h(ri,... ,r7v). (2) 


The ground-state wavefunction <h(ri,..., fat) is normalized to unity. We would like to ex¬ 
amine the many-body dynamics when the system is released from the trap 


H(ri, ..., r^v; t) = i 




dt 


T(ri,..., r^; f) = ..., r^). 


( 3 ) 


The initial condition is therefore T(ri,..., ftv; 0) = <h(ri,..., fat). 

To discuss the dynamics it is useful to make a coordinate transformation, from the lab¬ 
oratory frame ri,...,rAr to the center-of-mass {CM) and relative-motion {ret) coordinates 
Qi,..., Qat. Explicitly, we dehne (see, e.g.. Ref. [17]) 

k 




^/k{k + l) ^ 

, N 


^(Ffc+I-Fj), l<k<N-l, 




( 4 ) 
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The corresponding conjugated momenta are where, in particular, the momentum con¬ 
jugated to Qat equals Here, the respective Cartesian components are 

Qfc = {Qk,x, Qk,y, Qk,z) and ^ = (^g^, 0 ^, • With the help of the center-of-mass 

and relative-motion coordinates, the Hamiltonian after the trapping potential is switched 
off becomes separable and reads 




1 

Qvf-l, Qv) = —77- 


-HfeKQl) • • • ) Qvf—1 


(5) 


2dQj, 

where Hrei collects all other and only relative-motion terms. 

We can now proceed to the time evolution of the expectation values of interest. Let us 
look at the center-of-mass position and momentum operators in the x direction 


1 


N 


N 


( 6 ) 


PcM = y^^pj, 

i=i 

They are proportional, respectively, to the x component of Qat and its 


WcM - 

i=i 


1 d 


The center-of- 


where h,- = 

^ J I OXj 

conjugated momentum 7 ^, i.e., Xcm = -^Qn,x and Pcm = 
mass position and momentum operators satisfy the commutation relation [Xcm, Pcm] = i 
for any number of bosons N. 

We are interested in the evolution in time of the center-of-mass position-momentum 
uncertainty product, thus we need the respective variances 

,(t) = = Ay, ( 0 ) + 


A^- 

Xcm' 


{^\XcmPcm + PcmXcmI^) — 2 (< h | Ac - M |* h )(< h | PcA ^|$) 


1 A 2 

j. 1 Pcm 

N m 


( 0 ) 


• r 


A%Jt) = {^{t)\P^M\^{t)) - {^{t)\PcM\^m 




(7) 


where Eqs. (3,5,6) and the relations e~^'^^*XcMe~^^* = Xcm + and = 

Pcm have been used. We see that the variances at time t > 0 can be expressed by quantities 
at f = 0 only. This is a useful many-body relation. Note that the term linear in time can 
vanish on the following accounts: For the ground state in a general trap (<h|PcA^|<h) = 0; in 
a reflection-symmetric trap (<h|Af;A 4 '|‘h) = 0 ; and in an harmonic trap, see for the example 
below, {^\XcmPcm + PcmXcm\^) = 0. Nonetheless, we keep this term to account for the 
most general case. Finally, multiplying the center-of-mass variances in Fq. (7) we get 




Pcm 


Xcm ' 


PCA 


( 8 ) 


+ ( {^I^cmPcm + PcmXcmI^) — 2(<h| A(7 m Id*) (<h I Pcm I *h) 
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t + 
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for their uncertainty product. 

We are now in the comfortable position to examine a particular and solvable case of 
interest. The fact that one can provide an analytical result is attractive and instrumental in 
what follows, but by no means does it restrict the generality of the conclusions. We consider 
a system prepared in the ground state of the harmonic potential Vrir) = where u is 

the frequency. Since now ^(Qi,..., Qv-i, Qv) = • • • > Qa^-i); where 

^rei depends on the relative-motion coordinates only, the term linear in time of Eq. (8) 
vanishes and the other terms at t = 0 are readily evaluated and read 

(0) = (0) = —. (9) 

^CM ' ’ 2 uj N ^ 2 , ^ 

The hnal result for the center-of-mass position-momentum uncertainty product of an inter¬ 
acting system released out of the harmonic trap simplihes and is given by 

= ( 10 ) 
The result Eq. (10) holds for any number of particles N and general boson-boson interaction 
W{ri — r 2 ). In particular, it is the same behavior of a single free quantum particle initially- 
prepared in a Gaussian wavefunction which is left to propagate in time [1]. 

We would like to contrast the analytical, many-body result from what one would obtain 


from the Gross-Pitaevskii theory, which is expected to describe the exact out-of-equilibrium 
dynamics in the limit of inhnite number of particles [15]. When the Gross-Pitaevskii equa¬ 
tion, 

■ _ r \ d , \ / AT T\ f A* „I\ a a /I 1 \ 


i - 1 ) j dr’ -r’)(j)Gp{r-,t) (j)Gp{r-,t), (11) 

is employed to describe the system, i.e., for \hG'p(ri,... ,rN',t) = Y[f=i(pGp{'^j',t), the vari¬ 
ances of the center-of-mass position and momentum operators read 


{^Gp{t)\XcM\^Gp{t)) - {^Gp{t)\XcM\^Gpif)T = 

^ [(0Gp(^)|a)^|0Gp(^)) — (0Gp(O|a)|0Gp(^))^] , 

'^hcM,GpW ^ (^Gp{t)\PcM\^Gp{i)) - {^Gp{i)\PGM\^Gp{t)Y = 

= N [{(f)Gp{t)\p^\(l)Gp{t)) - (0gp(^)|p|?^igp(^))^] • (12) 


Gonsequently, the position-momentum uncertainty product at the mean-held level of theory. 


GP ^XcM,GpW^PcM,Gp(^) 

= [{(l>Gp{t)\x^\(i)Gp{t)) - (^Gp(t)|a:|0Gp(O)^] X [(0Gp(O|P^I^Gp(^)) “ (0Gp(O |p|^Gp(^))^] , 
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depends explicitly on the Gross-Pitaevskii orbital. 

Comparing Eqs. (10) and (13) implies that the time-dependent uncertainty product 
(t)A‘p^ (t) computed at the many-body level and at the Gross-Pitaevskii level cannot 
be the same. This conclusion is based on the shape of the time-dependent Gross-Pitaevskii 
orbital 4>Gp{r]t). Since the Gross-Pitaevskii orbital is distorted due to the particle-particle 
interaction from the shape of an expanding Gaussian (unless the interaction equals to 
zero), the position-momentum uncertainty product of any time-dependent non-Gaussian 
one-particle function is different than the minimal uncertainty product | (1 + ca^t^) of the 
exact result. This is a generalization of the result for the ground state [16] to the out-of- 
equilibrium dynamics of a freely expanding, initially-trapped BEG. 

In order to explore the size of the deviation of the Gross-Pitaevskii dynamics from the 
exact result we consider repulsive bosons in one spatial dimension released from a harmonic 
potential. The one-body Hamiltonian is — 2 ^ + the boson-boson interaction is con¬ 

tact, \qW{xi—X 2 ) = Ao(^(a;i—X 2 ). Since there is no analytical solution to the time-dependent 
Gross-Pitaevskii equation, Eq. (13), we solve it numerically. For the computation we have 
used 2000 sine discrete-variable-representation grid points in a box of size [—100,100]. Fig¬ 
ure 1 collects the results. In particular, the deviation of the mean-held result (13) from the 
analytical many-body relation (10) as a function of time is clearly visible. The stronger the 
interaction parameter A = Xq{N — 1) is, the larger the deviation is. This deviation can serve 
to quantify the amount of correlations in the system, and how this amount evolves in time 
during the dynamics. The many-body uncertainty product increases like with the same 
pace of growth for any interaction strength, see Eq. (10). The mean-held uncertainty prod¬ 
uct is seen to increase similarly in time, but the pace of growth depends on the interaction 
parameter A. This is because the initially-trapped atomic cloud is broader than that of the 
non-interacting system and it expands faster with increasing A [20]. All in all, the time- 
dependent, interaction-dressed shape of the atomic cloud enters directly the computation 
of the uncertainty product {t) when it is performed at the mean-held level. 

We remind that this behavior occurs even in the limit of inhnite number of particles, when 
the system becomes fully condensed. With this example, the analytical theoretical part of 
our work is concluded. We move to the uncertainty product in the many-body dynamics 
of an out-of-equilibrium trapped BEG, for which a numerical solution and investigation are 
necessary. 
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FIG. 1. (Color online) Time-dependent center-of-mass position-momentum uncertainty product 

(t)A‘^~ (t) of a BEC released from an harmonic trap. An illustrative numerical example with 

XcM PcM 

one-dimensional bosons. The one-body Hamiltonian is — and the interboson interaction 

is contact, AoIE(xi — X 2 ) = Ao5(a:i — X 2 ). Shown and compared as a function of time t are the 
Gross-Pitaevskii results for the interaction parameters A = \q{N — 1) = 1 (in red; dashed), A = 10 
(in green; dashed-dotted), and A = 100 (in blue; dashed-double-dotted) and the analytical, many- 
body result 1(1-1-1^) valid VA [Eq. (10)] (in black; full curve). The position-momentum uncertainty 
product computed at the Gross-Pitaevskii level differs from the analytical, many-body result. The 
difference increases upon increasing A, meaning that the pace of growth of the uncertainty at the 
mean-field level depends on the interaction parameter A (the y-axis is plotted in logarithmic scale). 
The many-body uncertainty product grows as and the mean-field uncertainty product is seen 
to grow in a similar manner in time (see Ref. [20] for the mean-field analysis of the expansion). 
The uncertainty product constitutes a macroscopic probe of the time-dependent correlations of a 
BEG, even when the system becomes 100% condensed in the limit of infinite number of particles. 
See the text for more details. The quantities shown are dimensionless. 







III. THE UNCERTAINTY PRODUCT IN THE DYNAMICS OF AN OUT-OF- 
EQUILIBRIUM TRAPPED BOSE-EINSTEIN CONDENSATE 

A. Theory 

The time-dependent Schrodinger equation of a trapped BEC, 

Hriri, ..., rjv)'h(ri,..., rjv; t) = (14) 

where the wavefunction \l/(ri,..., fat; t) is normalized to one, generally has no analyti¬ 
cal solution. Thus, a numerical solution of the out-of-equilibrium dynamics is usually a 
must. Furthermore and in the present context, there is, to the best of our knowledge, no 
mathematically-rigorous proof that the out-of-equilibrium many-body Schrodinger dynamics 
of a trapped BEC is described - in the limit of inhnite number of particles and at constant 
interaction parameter ~ by the time-dependent Gross-Pitaevskii equation. This is unlike 
the case of the ground state [13, 14] or that of an expanding, initially-trapped BEC [15]. 
Thus, the additional merit of a many-body numerical treatment here is to provide solid 
information on what can happen in this limit. 

It is useful to proceed by employing the reduced one-body and two-body density matrices 
of 4/(ri,..., r^r; t) [17, 21-24]. The reduced one-body density matrix is given by 

-—= j dr2.. ..., t)T(ri, ra, ...,rN,t) = 

= 5^^^ai(ri;t)a*(F'i;t). (15) 

j 

The quantities aj{r;t) are the so-called natural orbitals and Ujit) their respective occupa¬ 
tions which are time dependent and used to dehne the (possibly varying in time) degree 
of condensation in a system of interacting bosons [25]. The density of the system is the 
diagonal of the reduced one-body density matrix, p(r; f) = p‘^^^(r, r; t). 

It is convenient to express in what follows quantities using the time-dependent natural 
orbitals aj{r-, t). The diagonal part of the reduced two-body density matrix is given by 

p(2 )(ri,r2,ri,r2;t) f 

-- = J dr3...(irjvT (fi, r2 ,..., fjv; t)T(ri, r2 ,..., r^v; t) = 

= mT-\) »*j{^i,t)a*p{r2-,t)ak{ri,t)aq{r2-,t), (16) 

jpkq 
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where the matrix elements are pjpkqif) = {^{t)Wjbl,bkbg\'^{t)), and the creation and annihi¬ 
lation operators are associated with the time-dependent natural orbitals. 

We are interested in extracting the uncertainty product of operators for the time- 
dependent state of N interacting bosons in a trap described by the wavefunction 
'l/(ri,..., r^r; f). We begin simply with the operator 

N 

i = ^a(rj) (17) 

i=i 

of the many-particle system, where a(r) is an hermitian operator. A straightforward calcula¬ 
tion gives the average per particle of A in the state ;^(\['(t)|A|d'(t)) = J dr^^^^a(r), 

which is seen to be directly related to the density of the system [26]. 

To compute the variance we also need the expectation value of the square of A, 

N 

A^ = ^cA{rj)+ '^ 2 a{rj)d{rk), (18) 

j=l j<k 

which is comprised of one-body and two-body operators. Expressed in terms of the above 
quantities, the time-dependent variance per particle of the operator A of the many-particle 
system reads 


1 

N 



i ((>p(i)q"i'i'(i)) ■ 

- (*(()|i|*(i)>") 

f J PAA) 2f N 

J dr ^ a (r) 

1-1 

u 


N 


1 2 


dr|Q;i(r; t)\‘^a{ 


(Af-l) 



dr 


N 



+ 




Pjpkqif) 


N 


dra*,{r] t)(^kA') t)A 


dra*{r] t)aq{r] t)a{ 


(19) 


jpkq^llll 

where we reordered the terms as explained below. The hrst two describe the variance of a(r) 
computed using the density per particle The next two terms essentially cancel each 

other for condensed systems. They exactly cancel each other within Gross-Pitaevskii theory, 
since then piiii(t) = N{N — 1) and \ai{r;t)\‘^ = |^Gp(i‘;t)P = The last sum can 

contribute only when a theoretical description beyond Gross-Pitaevskii is employed. This 
allows the variance A^(f) to pick-up information from the other terms of the reduced two- 
body density matrix, PjpkqiA, jpkq A HH, that normally do not contribute. Furthermore, 
the explicit dependence of the last term in Eq. (19) on integrals involving the time-dependent 
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natural orbitals {ak{r-,t)} means that the contribution of many-body effects depends also 
on the instantaneous shape of the natural orbitals. We will return to and examine this issue 
below. Finally, in the absence of interboson interaction, the time-dependent variance per 
particle boils down to that of a single particle. 

It is useful for the sake of subsequent analysis to look at the difference in the variance of 
a system when it is described at the many-body and Gross-Pitaevskii levels: 


^ = ^Icorrelationsit), 


N 


N A,GP^ 



f *h(r)l 

J N ^ ^ 

J N ^ 


( 20 ) 


We henceforth call this difference the correlations term, A\ . The correlations 

term vanishes identically for non-interacting bosons, when all the bosons occupy one and the 
same time-dependent orbital. For interacting bosons, if natural orbitals other than the first 
natural orbital (the so-called condensed mode) are at all occupied, A\ does not 

vanish, see also Eq. (19). To remind, the depletion of the BEG is given by the occupation of 
all but the first natural orbital divided by N. Eor condensed systems the depletion is very 
small. So, the question is whether a tiny occupation of the higher natural orbitals can make 
a difference? We point out that in an out-of-equilibrium dynamics this occupation may vary 
in time and still remain tiny. We already demonstrated in Sec. II that the answer is positive 
for the freely expanding BEG. We shall see below that the answer remains positive for the 
dynamics of a trapped BEG as well. 


Next, consider two many-particle operators, Eq. (17), A = B = 


'n!j=i Kg)’ their respective time-dependent variances per particle, ^A^(t) and ^A|(f). 


The uncertainty product of these two operators satishes the inequality [26] 

|2 


AT ' N ^ 


N 


*d^[a(r),i)(r)l 


( 21 ) 


where [A, B] = A[a(r), b(r)] and [a(r), 6(r)] is the usual commutator of the hermitian oper¬ 
ators a(r) and b{r). Using Eq. (20) we can express the difference in the uncertainty product 
when it is described at the many-body and Gross-Pitaevskii levels as follows: 


= A^ 

^decorrelations 




+A 


d,correlations d)^i,,correlations^^^ ’ 


( 22 ) 
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We see that the difference in the uncertainty product Eq. (22) consists of contributions from 
the correlations term in the variance of A times the mean-held variance of B, from the vice 
versa expression, and from the product of the correlations terms in the variances of A and 
B. We will return to this analysis in the numerical study below. 


B. Case study: A bosonic Josephson junction 

One of the most studied systems in the context of trapped BECs is the bosonic Joseph¬ 
son junction [27-39]. We therefore would like to implement the above ideas and investigate 
numerically the time-dependent center-of-mass position-momentum uncertainty product of 
a BEG in a double-well potential. We need a suitable and proved many-body tool to arrive 
at detailed conclusions. Such a many-body tool is the multiconhgurational time-dependent 
Hartree for bosons (MCTDHB) method, which has been well documented [40-43], bench- 
marked [44], and extensively applied [45-59] in the literature. Particularly, the numerically- 
exact quantum dynamics of the one-dimensional bosonic-Josephson-junction system has 
been reported in Ref. [36]. 

Consider a symmetric double-well potential Vt{x) formed by fusing together the two 
“left” and “right” harmonic potentials |(x ± 2)^ (see the appendix for details). We begin 
by preparing at t < 0 a BEG made of N interacting bosons in the ground state of the left 
harmonic trap. At t > 0, the system is let to evolve in time in the double-well trap Vt{x). 
The bosons-boson interaction is contact, \qW{xi — X 2 ) = AoJ(xi — X 2 ). The interaction 
parameter is A = \q{N — 1) = 0.01. 

The time-dependent survival probability in the left well, Piit) = , of the BEG 

in the one-dimensional bosonic Josephson junction is registered in Fig. 2. We compare the 
many-body time evolution of A = 1 000, 10 000, and 100 000 bosons using M = 2 time- 
dependent orbitals, and the Gross-Pitaevskii, mean-held result (M = 1). It is found that all 
four curves, the three many-body calculations and the mean-held one, he atop each other. 
We observe full tunneling of the density back and forth between the left and right wells. 
The slight deviation from the maximal (minimal) value of 100% (0%) survival probability, 
and the small beating of the maxima of Piit), are due to the hnite depth of the double-well 
and the preparation of the initial condition in the left harmonic potential. The apparent 
discontinuities are in fact smooth, see the inset in Fig. 2. Important to our investigations 
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^ ^Rabi 

FIG. 2. (Color online) (a) Time-dependent survival probability in the left well, Piit), of a BEC in a 
one-dimensional bosonic Josephson junction. The interboson interaction is contact, \qW{xi—X 2 ) = 
Ao(5(xi — X 2 )- The interaction parameter is A = \q{N — 1) = 0.01. Shown are the many-body 
results for N = 1 000, 10 000, and 100 000 bosons (in color; full curves) using M = 2 time-dependent 
orbitals, and the Gross-Pitaevskii, mean-field result (in black; dashed curve). All four curves, the 
three many-body results and the mean-field one, lie on top of each other. Full tunneling of the 
density back and forth between the left and right wells in seen. The inset shows the smoothness 
of the survival probability at its extrema. The panels to the right depict snapshots of the time- 
dependent density (with the double-well potential in the background) at three different times, when 
the density is (b) either localized in the left well, (c) delocalized over both wells, (d) or localized 
in the right well. The many-body (in color; full curves) and mean-field (in black; dashed curve) 
results lie on top of each other. That the many-body and mean-field results for the time-dependent 
density and survival probability coincide indicates that the system is condensed (also see Fig. 3b). 
See the text for more details. The quantities shown are dimensionless. 

is that the many-body and mean-held results for the time-dependent density and survival 
probability coincide. This indicates that the system is fully condensed (also see Fig. 3b in 
this respect). Furthermore, the interaction is so weak that known many-body effects in the 
bosonic Josephson junction, such as the collapse of the density oscillations [27, 38], have not 
yet developed. 
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Fig. 3a plots the time-dependent center-of-mass position-momentum uncertainty product 
(t)A‘p^ (t) in the bosonic Josephson junction [60]. The time-dependent uncertainty 
product computed at the many-body level is seen to grow in time in an oscillatory manner. 
In contrast to the many-body result, the Gross-Pitaevskii uncertainty product only oscillates 
and thus deviates from the many-body results as time progresses. This represents a major 
difference between the mean-held and many-body dynamics of the BEG. Note that the results 
for N = 1000, 10 000, and 100 000 bosons he on top of each other and suggest convergence 
of the many-body uncertainty product with increasing number of bosons and at constant 
interaction parameter all the way to the inhnite particle limit. 

Along with the oscillatory growth of the uncertainty product, the depletion of the (hnite) 
BEG starts to increase, see Fig. 3b. We observe that it is enough to have less than one 
particle outside the condensed mode in order for the many-body uncertainty product to 
macroscopically differ from the mean-held one. Importantly, one can also see that the 
time-dependent depletion decreases with increasing number of bosons N and at constant 
interaction parameter A. This suggests that the trapped time-dependent BEG remains 100% 
condensed in the limit of inhnite number of particles. Interestingly, the curves in Fig. 3b 
appear to be shifted vertically from each other, suggesting that the total time-dependent 
number of particles residing in the higher natural orbitals does not depend on N. The latter 
time-dependent observation is analogous to the static situation for the ground state, see 
Fig. Ic in Ref. [16]. 

Let us analyze the dynamics of the uncertainty product more closely. There are several 
features seen in Fig. 3a that we address: The period of oscillations, the minimal and maximal 
values, the hne structure atop the curves and, prominently, the growth pace itself. Since the 
BEG is fully condensed, it is instrumental to hrst discuss the mean-held dynamics and then 
the many-body one. Because the repulsion is very weak, the period of oscillations of Piit) is 
essentially that of the non-interacting system tRabi, see Fig. 2. In each oscillations’ cycle the 
BEG is localized once in the left and once in the right well. Hence, the period of oscillations 
of the uncertainty product is half the period of the survival probability, i.e., \tRaU- 

The minimal and maximal values of the center-of-mass uncertainty product can for the 
Gross-Pitaevskii dynamics be estimated as follows, making use that the repulsion is very 
weak. When the BEG is localized in either the left or right well, the minimal uncertainty 
product is essentially that of the initial-condition mean-held ground-state in the harmonic 
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FIG. 3. (Color online) (a) Time-dependent center-of-mass position-momentum uncertainty product 
(t) of a BEG in the one-dimensional bosonic Josephson junction of Fig. 2. Shown 

^CM PcM 

are many-body results using M = 2 time-dependent orbitals. The interaction parameter is A = 
Xo{N — 1) = 0.01. The system is weakly interacting. The time-dependent uncertainty product 
computed at the many-body level is seen to grow in time in an oscillatory manner. The period of 
oscillations is essentially \tRabi- For comparison and later analysis, the Gross-Pitaevskii dynamics 
only oscillates and thus seen to deviate from the many-body results as time progresses. The 
results for N = 1000 (in red), 10000 (in green), and 100000 (in light blue) bosons lie on top 
of each other and suggest convergence of the many-body uncertainty product with increasing 
number of bosons and at constant interaction parameter in the infinite particle limit, (b) The 
time-dependent depletion decreases with N and at constant A (note the small values on the y axis 
which is in logarithmic scale). Side by side, the curves appear to be shifted vertically from each 
other, suggesting that the total time-dependent number of particles residing in the higher natural 
orbitals is constant. See the text for more details. The quantities shown are dimensionless. 
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potential, i.e., r^r, ~ w ‘ 4 = i- the other hand, when the BEC is 

delocalized over the two wells and assuming their inter-well distance to be large (the inter¬ 
well distance is equal to 4), the center-of-mass position variance is (I + ~ 

^ and the momentum variance is still about Ap ^ y. Consequently, the maximal 
uncertainty product is A^^ gp^Pc gp ^ W ' T ^ i- These extreme values of the 
uncertainty product nicely match the mean-held numerical results in Fig. 3. The small 
deviations from these analytical values, as well as the hue structure atop the maximal values 
of the uncertainty product, are due to the hnite depth and size of the double well, which lead 
to a small mixture between the wells. Consequently, the initial condition prepared in the 
left harmonic well slightly penetrates the right well, and the sudden change at f = 0 from a 
single to a double well along with the weak interaction mildly couple the initial wavepacket 
to higher than the lowest two modes of the non-interacting double well. 

We are now in the position to discuss the many-body uncertainty. At the very beginning 
of the density oscillations, the mean-held and many-body uncertainty products he atop each 
other, see Fig. 3a. Then they start to deviate from each other in a unique manner. The 
maxima of the many-body uncertainty product in the double well starts to grow more and 
more in time (and the hue structure atop the curves is magnihed), however the minima in 
the uncertainty product essentially coincide. Also the period of oscillations are the same, see 
the respective survival probabilities in Fig. 2. What is the reason for this unique difference? 

Since the many-body and mean-held densities coincide and the system is condensed, see 
Figs. 2 and 3b, respectively, the correlations term of the variance equals essentially only 
the last sum in (19), i.e., ^ T,jpkq^uu^^^ [J dxa*{x;t)ak{x;t)a{x)] x 

[/ dxa*{x;t)ag{x;t)a{x)~\. Even the slightest occupation of the non-condensed modes, 
see Fig. 3b, is enough to generate a sizable value of the correlations term and 
this macroscopically impacts the many-body variance, see Fig. 3a. 

We are left now to analyze the individual contributions of the position A\ 
and momentum correlations terms to the dynamics of the uncertainty prod¬ 

uct in the bosonic-Josephson-junction system. A key point to note is that the density 
oscillates between the left and right wells back and forth. Furthermore, the interaction 
is so weak that two time-dependent natural orbitals are sufficient to converge the dynam¬ 
ics of the uncertainty product, please see Fig. 5 in the appendix. When the density is 
localized in, say, the left well, both natural orbitals are localized in this well, and the 
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various contributions of the overlap integrals between the natural orbitals times the re¬ 
spective elements of the reduced two-body density matrix to the position and momentum 
correlations terms are small. On the other hand, when the density is spread over the two 
wells, both natural orbitals, the condensed ai{x]t) and excited a 2 {x]t) modes, are de¬ 
localized over the junction and mimic gerade and ungerade shapes. Then, only the ele¬ 
ment [/ dxal {x;t)a 2 {x;t)x] [J dxal{x;t)a 2 {x;t)x] has a dominant contribution to 

the position correlations term. Finally, examining the momentum correlations term, the 
only possible dominant element, [/ dxal{x;t)^-^a2{x;t)~\ [J dxal{x;t)^-^a2{x;t)~\, 

is, however, small as well. This is because the delocalized natural orbitals look like linear 
combinations of Gaussian-like functions, and taking their spatial derivatives leads to func¬ 
tions having nodes in each well, and thus small overlaps with the natural orbitals themselves. 

Combining all the above, when the density is spread over the two wells, the elements of the 
reduced two-body density matrix times the overlap integrals containing the natural orbitals 
give a dominant many-body contribution to the uncertainty product. Furthermore, analysis 
of the instantaneous shape of the time-dependent natural orbitals traces this contribution to 
^Icorreiationsit) rather than Interestingly, the many-body position variance 

surpasses the trap’s size, which sets in the maximal value of mean-held position variance. 
This might be perceived as surprising, because we are dealing with a condensed system 
conhned in a trap. But it is not, being a genuine many-body effect of the interacting 
particles. Summing up, the main contribution to the difference between the many-body and 
mean-held uncertainty product, see Eq. (22), in the bosonic Josephson junction dynamics 
comes from the product of the time-dependent position correlations term and the mildly 
changing mean-held variance of the momentum operator itself. 

Finally, we set to examine the growth of the many-body uncertainty product in time. 
In Fig. 4 we have computed the many-body dynamics of = 1000 bosons for the three 
interaction parameters A = 0.005, 0.01, and 0.02. The larger the interaction parameter 
is, the faster the many-body uncertainty product grows (in an oscillatory manner). We 
concentrate on the maxima values of the many-body uncertainty, since at the minima the 
many-body and mean-held results practically coincide. A least-squares ht to a polynom of 
the form a(A) -|- b{A) ■ was made to the maxima of each of the three curves in Fig. 4. 
We can hrst see that the square growth in time nicely matches the maxima values of the 
uncertainty. We can numerically conclude that, to leading order in time, the (oscillating) 
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FIG. 4. (Color online) Time-dependent center-of-mass position-momentum uncertainty product 
(t) of a BEC in the one-dimensional bosonic Josephson junction of Fig. 2: Analysis 

■^CM ^CM 

of its growth in time. Results for N = 1 000 bosons for the interaction parameters A = Ao(A^ — 1) = 
0.005, 0.01, and 0.02 are depicted. Shown are many-body results using M = 2 time-dependent 
orbitals. The black curves are least-squares hts of the maxima of the uncertainty product to 
a polynom of the form a(A) -|- 6 (A) • . The numerical results suggest that the envelop of the 

time-dependent uncertainty product grows as A^t^, to leading orders in time and the interaction 
parameter. See the text for more details. The quantities shown are dimensionless. 


many-body uncertainty product of the BEC grows in the bosonic Josephson junction like 

We can learn more from examining the interaction-dependent polynomial’s coefficients 
a(A) and 6(A). The respective values are a(0.005) = 2.186 and 6(0.005) = 1.749 • 10“®; 
a(O.Ol) = 2.153 and 6(0.01) = 6.953 • 10“®; and a(0.02) = 2.019 and 6(0.02) = 2.782 • 
10“®. The baseline values are nearly the same, and account for the many-body value of 
the uncertainty of the BEC, when it is momentarily delocalized between the two wells 
(extrapolated to the beginning of the propagation). These values are close to but slightly 
smaller than the analytical value of | discussed above within the analysis of the mean- 
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field uncertainty product. Furthermore, we find the following ratios for different interaction 
parameters: = 3.976 ~ 4 and = 4.002 Ri 4. We can numerically conclude that, 

to leading order in the interaction parameter A, the (oscillating) many-body uncertainty 
product of the BEC in the bosonic Josephson junction grows like A^. Consequently and 
all together, the many-body uncertainty prodnct in the double well grows to leading orders 
nnmerically like A^f^. Finding an analytical explanation for this nnmerically-fonnd growth 
law of the nncertainty prodnct would be interesting bnt goes beyond the scope of the present 
work. We speculate that Ref. [39] could be generalized and instrnmental here. Of course, 
hnite BECs in the bosonic Josephson jnnction are interesting for their own right. Then, as 
time goes by or when the interaction strength is increased, known many-body effects, like 
the collapse of the density oscillations [27, 38], eventnally set in. It wonld be interesting to 
explore the nncertainty prodnct when the system is no longer condensed. 


IV. CONCLUDING REMARKS 

We have discussed in the present work, analytically and nnmerically, the time-dependent 
nncertainty prodnct of an ont-of-eqnilibrinm BEC. Althongh the systems considered are 
fnlly condensed, the uncertainty prodnct at the many-body level differs from that predicted 
at the Gross-Pitaevskii, mean-held level. As the system evolves in time, this difference can 
become snbstantial. In this capacity, onr work demonstrates that one must use a many-body 
propagation theory to describe the dynamics of an ont-of-eqnilibrinm BEC, even in the limit 
of inhnite nnmber of particles when the system becomes 100% condensed. 

To measnre this time-dependent many-body effect, the position and momentnm of all 
particles in the clonds are in principle needed. This is in contrast to the commonly reported 
particle density. Experimental techniqnes which conld achieve the necessary resolntion are 
constantly being improved [63-67]. 
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FIG. 5. (Color online) Time-dependent center-of-mass position-momentum uncertainty product 
{t)A‘p (t) of a BEG in the one-dimensional bosonic Josephson junction of Fig. 2: Conver¬ 
gence of the out-of-equilibrium many-body dynamics. Results for N = 100 bosons for an increasing 
number M of the time-dependent orbitals are shown. The curves lie on top of each other, and the 
uncertainty product seen to converge with increasing M. The quantities shown are dimensionless. 


APPENDIX: DETAILS OF THE TIME-DEPENDENT NUMERICAL COMPU¬ 
TATIONS 

The multiconfigurational time-dependent Hartree for bosons (MCTDHB) method [40-43] 
utilizes M time-dependent orbitals, which are determined according to the Dirac-Frenkel, 
time-dependent variational principle, to solve the time-dependent many-boson Schrodinger 
equation Hriri ,..., r 7 v)tl/(ri,..., r^v; t) = As the method is variational, con¬ 

vergences in the limit M —)■ oo orbitals to the solution of the time-dependent many-boson 
Schrodinger equation is achieved [44]. On the other end, i.e., for M = 1, the MCTDHB 
method boils down to the time-dependent Gross-Pitaevskii, mean-held theory. The employ¬ 
ment of M time-dependent orbitals allows one to obtain accurate numerical results with 
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substantially less numerical resources than a corresponding time-dependent computation 
employing M fixed orbitals would need. To obtain the initial conditions for the time prop¬ 
agations, which here are the respective many-body ground states in the harmonic trap, we 
propagate the MCTDHB equations-of-motion in imaginary time [44, 61]. The MCTDHB 
method has recently been benchmarked and its accuracy demonstrated [44]. We use the 
numerical implementation in the software package [62]. 

In the present work the MCTDHB method is employed to compute the time-dependent 
uncertainty product in the out-of-equilibrium dynamics of a one-dimensional bosonic Joseph- 
son junction. The many-body Hamiltonian is represented by 201 harmonic-oscillator 
discrete-variable-representation grid points in a box of size [-10,10]. The symmetric double¬ 
well potential is obtained by merging the two harmonic potentials |(x ± 2)^, whose minima 
are located, respectively, to the left and right of the origin at a; = =f 2, with a square- 
polynomial-barrier in the central region |x| < giving rise to Vt{x) = + 2)^,a; < 

— |(1 — a;^), |x| < |(a: — 2)^, a; > |}. A convenient unit of time in the double well is 

the tunneling time, fRabi = ^ = 132.498, where Ai? is the energy difference between the 
ground state and the first excited state of a single particle in the double-well potential. 

As a concrete example and without loss of generality, convergence with increasing M 
of the time-dependent many-particle position-momentum uncertainty product of an out- 
of-equilibrium BEC held in the one-dimensional bosonic Josephson junction of Fig. 2 is 
demonstrated in Fig. 5. Recall that the density and survival probability converge already 
and the Gross-Pitaevskii level (M = 1), see Fig. 2. Fig. 5 shows that the dynamics of the 
many-particle position-momentum uncertainty product converges at the M = 2 level of the 
MCTDHB theory. 


[1] C. Cohen-Tannoudji, B. Diu, and F. Laloe, Quantum Mechanics, Vol. 1 (Wiley, New York, 
1977). 

[2] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Observation 
of Bose-Einstein Condensation in a Dilute Atomic Vapor, Science 269 , 198 (1995). 

[3] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Evidence of Bose-Einstein Con¬ 
densation in an Atomic Gas with Attractive Interactions, Phys. Rev. Lett. 75 , 1687 (1995). 


21 



[4] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and 
W. Ketterle, Bose-Einstein Condensation in a Gas of Sodium Atoms, Phys. Rev. Lett. 75, 
3969 (1995). 

[5] E. A. Cornell and C. E. Wieman, Nobel Lecture: Bose-Einstein condensation in a dilute gas, 
the first 70 years and some recent experiments, Rev. Mod. Phys. 74 , 875 (2002). 

[6] W. Ketterle, Nobel lecture; When atoms behave as waves: Bose-Einstein condensation and 
the atom laser. Rev. Mod. Phys. 74 , 1131 (2002). 

[7] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of Bose-Einstein condensa¬ 
tion in trapped gases, Rev. Mod. Phys. 71 , 463 (1999). 

[8] A. J. Leggett, Bose-Einstein condensation in the alkali gases: Some fundamental concepts, 
Rev. Mod. Phys. 73 , 307 (2001). 

[9] 1. Bloch, J. Dalibard, and W. Zwerger, Many-body physics with ultracold gases Rev. Mod. 
Phys. 80 , 885 (2008). 

[10] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford University Press, Oxford, 
2003). 

[11] A. J. Leggett, Quantum Liquids: Bose condensation and Cooper pairing in condensed matter 
systems (Oxford University Press, Oxford, 2006). 

[12] G. J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases, 2nd ed. (Gambridge 
University Press, Cambridge, England, 2008). 

[13] E. H. Lieb, R. Seiringer, and J. Yngvason, Bosons in a trap: A rigorous derivation of the 
Gross-Pitaevskii energy functional, Phys. Rev. A 61 , 043602 (2000). 

[14] E. H. Lieb and R. Seiringer, Proof of Bose-Einstein Gondensation for Dilute Trapped Gases, 
Phys. Rev. Lett. 88, 170409 (2002). 

[15] L. Erdos, B. Schlein, and H.-T. Yau, Rigorous Derivation of the Gross-Pitaevskii Equation, 
Phys. Rev. Lett. 98, 040404 (2007). 

[16] S. Klaiman and O. E. Alon, Variance as a sensitive probe of correlations, Phys. Rev. A 91 , 
063613 (2015). 

[17] L. Cohen and C. Lee, Exact reduced density matrices for a model problem, J. Math. Phys. 
26 , 3105 (1985). 

[18] Y. Castin and R. Dum, Bose-Einstein Condensates in Time Dependent Traps, Phys. Rev. 
Lett. 77 , 5315 (1996). 


22 



[19] L. Plaja and L. Santos, Expansion of a Bose-Einstein condensate in an atomic waveguide, 
Phys. Rev. A 65 , 035602 (2002). 

[20] V. A. Brazhnyi, A. M. Kamchatnov, and V. V. Konotop, Hydrodynamic flow of expanding 
Bose-Einstein condensates, Phys. Rev. A 68, 035603 (2003). 

[21] P.-O. Lowdin, Quantum Theory of Many-Particle Systems. I. Physical Interpretations by 
Means of Density Matrices, Natural Spin-Orbitals, and Convergence Problems in the Method 
of Configurational Interaction, Phys. Rev. 97 , 1474 (1955). 

[22] A. J. Coleman and V. I. Yukalov, Reduced Density Matrices: Coulson’s Challenge, Lectures 
Notes in Chemistry Vol. 72 (Springer, Berlin, 2000). 

[23] Reduced-Density-Matrix Mechanics: with Application to Many-Electron Atoms and Molecules, 
edited by D. A. Mazziotti, Advances in Chemical Physics Vol. 134 (Wiley, New York, 2007). 

[24] K. Sakmann, A. I. Streltsov, O. E. Alon, and L. S. Cederbaum, Reduced density matrices and 
coherence of trapped interacting bosons, Phys. Rev. A 78 , 023615 (2008). 

[25] O. Penrose and L. Onsager, Bose-Einstein Condensation and Liquid Helium, Phys. Rev. 104 , 
576 (1956). 

[26] This is in case of a local operator in coordinate space r. In the respective case in momentum 
space p, ^(T(t)l A|'I'(t)) = f dp ^^^’*^ a(p). Por a ‘mixed’ operator, like the angular momentum 
operator, one generally requires its matrix elements with respect to the time-dependent natural 
orbitals. 

[27] G. J. Milburn, J. Corney, E. M. Wright, and D. E. Walls, Quantum dynamics of an atomic 
Bose-Einstein condensate in a double-well potential, Phys. Rev. A 55, 4318 (1997). 

[28] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, Quantum Coherent Atomic Tunneling 
between Two Trapped Bose-Einstein Condensates, Phys. Rev. Lett. 79 , 4950 (1997). 

[29] S. Raghavan, A. Smerzi, and V. M. Kenkre, Transitions in coherent oscillations between two 
trapped Bose-Einstein condensates, Phys. Rev. A 60 , R1787 (1999). 

[30] C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, and M. A. Kasevich, Squeezed States 
in a Bose-Einstein Condensate, Science 291 , 2386 (2001). 

[31] A. Vardi and J. R. Anglin, Bose-Einstein Condensates beyond Mean Field Theory; Quantum 
Backreaction as Decoherence, Phys. Rev. Lett. 86, 568 (2001). 

[32] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Gristiani, and M. K. Oberthaler, Direct 
Observation of Tunneling and Nonlinear Self-Trapping in a Single Bosonic Josephson Junction, 


23 



Phys. Rev. Lett. 95 , 010402 (2005). 

[33] T. Schumm, S. HofFerberth, L. M. Andersson, S. Wildermuth, S. Groth, I. Bar-Joseph, J. 
Schmiedmayer, and P. Kruger, Matter-wave interferometry in a double well on an atom chip, 
Nat. Phys. 1 , 57 (2005). 

[34] S. Levy, E. Lahoud, 1. Shomroni, and J. Steinhauer, The a.c. and d.c. Josephson effects in a 
BoseEinstein condensate. Nature (London) 449 , 579 (2007). 

[35] M. Trujillo-Martinez, A. Posazhennikova, and J. Kroha, Nonequilibrium Josephson Oscilla¬ 
tions in Bose-Einstein Condensates without Dissipation, Phys. Rev. Lett. 103 , 105302 (2009). 

[36] K. Sakmann, A. 1. Streltsov, O. E. Alon, and L. S. Cederbaum, Exact Quantum Dynamics of 
a Bosonic Josephson Junction, Phys. Rev. Lett. 103 , 220601 (2009). 

[37] T. Zibold, E. Nicklas, C. Gross, and M. K. Oberthaler, Classical Bifurcation at the Transition 
from Rabi to Josephson Dynamics, Phys. Rev. Lett. 105 , 204101 (2010). 

[38] K. Sakmann, A. 1. Streltsov, O. E. Alon, and L. S. Cederbaum, Universality of fragmentation 
in the Schrodinger dynamics of bosonic Josephson junctions, Phys. Rev. A 89 , 023602 (2014). 

[39] H. Veksler and S. Fishman, Semiclassical analysis of Bose-Hubbard dynamics. New J. Phys. 
17 , 053030 (2015). 

[40] A. 1. Streltsov, O. E. Alon, and L. S. Cederbaum, Role of Excited States in the Splitting of 
a Trapped Interacting Bose-Einstein Condensate by a Time-Dependent Barrier, Phys. Rev. 
Lett. 99 , 030402 (2007). 

[41] O. E. Alon, A. 1. Streltsov, and L. S. Cederbaum, Multiconfigurational time-dependent Hartree 
method for bosons: Many-body dynamics of bosonic systems, Phys. Rev. A 77 , 033613 (2008). 

[42] Multidimensional Quantum Dynamics: MCTDH Theory and Applications^ edited by H.-D. 
Meyer, F. Gatti, and G. A. Worth (Wiley-VCH, Weinheim, 2009). 

[43] Quantum Gases: Finite Temperature and Non-Equilibrium Dynamics, edited by N. P. 
Proukakis, S. A. Gardiner, M. J. Davis, and M. H. Szymanska, Cold Atoms Series Vol. 1 
(Imperial College Press, London, 2013). 

[44] A. U. J. Lode, K. Sakmann, O. E. Alon, L. S. Cederbaum, and A. 1. Streltsov, Numerically 
exact quantum dynamics of bosons with time-dependent interactions of harmonic type, Phys. 
Rev. A 86, 063606 (2012). 

[45] J. Grond, J. Schmiedmayer, and U. Hohenester, Optimizing number squeezing when splitting 
a mesoscopic condensate, Phys. Rev. A 79 , 021603(R) (2009). 


24 



[46] J. Grond, T. Betz, U. Hohenester, N. J. Mauser, J. Schmiedmayer, and T. Schumm, The 
Shapiro effect in atomchip-based bosonic Josephson junctions. New J. Phys. 13 , 065026 (2011). 

[47] M. Heimsoth, D. Hochstuhl, C. E. Crefheld, L. D. Carr, and F. Sols, Effective Josephson 
dynamics in resonantly driven Bose-Einstein condensates. New J. Phys. 15 , 103006 (2013). 

[48] R. Schmitz, S. Kronke, L. Cao, and P. Schmelcher, Quantum breathing dynamics of ultracold 
bosons in one-dimensional harmonic traps: Unraveling the pathway from few- to many-body 
systems, Phys. Rev. A 88 043601 (2013). 

[49] A. U. J. Lode, S. Klaiman, O. E. Alon, A. I. Streltsov, and L. S. Cederbaum, Controlling the 
velocities and the number of emitted particles in the tunneling to open space dynamics, Phys. 
Rev. A 89 , 053620 (2014). 

[50] O. I. Streltsova, O. E. Alon, L. S. Cederbaum, and A. I. Streltsov, Generic regimes of quantum 
many-body dynamics of trapped bosonic systems with strong repulsive interactions, Phys. Rev. 
A 89, 061602(R) (2014). 

[51] S. Klaiman, A. U. J. Lode, A. 1. Streltsov, L. S. Cederbaum, and O. E. Alon, Breaking the 
resilience of a two-dimensional Bose-Einstein condensate to fragmentation, Phys. Rev. A 90 , 
043620 (2014). 

[52] S. 1. Mistakidis, L. Cao, and P. Schmelcher, Negative-quench-induced excitation dynamics for 
ultracold bosons in one-dimensional lattices, Phys. Rev. A 91 033611 (2015). 

[53] S. Kronke and P. Schmelcher, Many-body processes in black and gray matter-wave solitons, 
Phys. Rev. A 91 053614 (2015). 

[54] U. R. Fischer, A. U. J. Lode, and B. Chatterjee, Condensate fragmentation as a sensitive 
measure of the quantum many-body behavior of bosons with long-range interactions, Phys. 
Rev. A 91 , 063621 (2015). 

[55] S. Kronke and P. Schmelcher, Two-body correlations and natural-orbital tomography in ul¬ 
tracold bosonic systems of definite parity, Phys. Rev. A 92 , 023631 (2015). 

[56] S. E. Weiner, M. C. Tsatsos, L. S. Cederbaum, and A. U. J. Lode, Angular momentum in in¬ 
teracting many-body systems hides in phantom vortices, arXiv;1409.7670v2 [cond-mat.quant- 
gas]. 

[57] 1. Brouzos, A. 1. Streltsov, A. Negretti, R. S. Said, T. Caneva, S. Montangero, and T. Calarco, 
Quantum Speed Limit and Optimal Control of Many-Boson Dynamics, arXiv;1412.6142v2 
[quant-ph]. 


25 



[58] K. Sakmann and M. Kasevich, Single shot simulations of dynamic quantum many-body sys¬ 
tems, arXiv:1501.03224vl [cond-mat.quant-gas]. 

[59] S. I. Mistakidis, T. Wulf, A. Negretti, and P. Schmelcher, Resonant quantum dynamics of few 
ultracold bosons in periodically driven finite lattices, arXiv:1506.01261vl [cond-mat.quant- 
gas]. 

[60] Note that Xqm = w Pcm = P implies that (^) = This 

means that the uncertainty product can be analyzed equivalently in terms of the variances 

it) and A|, it) or the variances per particle pA^Jt) and iA^(t). 

XcM rCM X ' iV ' 

[61] A. I. Streltsov, O. E. Alon, and L. S. Cederbaum, General variational many-body theory with 
complete self-consistency for trapped bosonic systems, Phys. Rev. A 73 , 063626 (2006). 

[62] A. I. Streltsov, L. S. Cederbaum, O. E. Alon, K. Sakmann, A. U. J. Lode, J. Grond, O. I. 
Streltsova, and S. Klaiman, The Multiconfigurational Time-Dependent Hartree for Bosons 
Package, version 3.x, http://mctdhb.org, Heidelberg (2006-Present). 

[63] W. S. Bakr, J. I. Gillen, A. Peng, S. Polling, and M. Greiner, A quantum gas microscope for 
detecting single atoms in a Hubbard-regime optical lattice. Nature (London) 462 , 74 (2009). 

[64] J. P. Sherson, C. Weitenberg, M. Endres, M. Cheneau, 1. Bloch, and S. Kuhr, Single-atom- 
resolved fluorescence imaging of an atomic Mott insulator, Nature (London) 467, 68 (2010). 

[65] W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. 1. Gillen, S. Foiling, L. Pollet, and 
M. Greiner, Probing the Superfluid-to~Mott Insulator Transition at the Single-Atom Level, 
Science 329 , 547 (2010). 

[66] M. Endres, M. Cheneau, T. Fukuhara, C. Weitenberg, P. SchauB, C. Gross, L. Mazza, M. 
C. Bahuls, L. Pollet, 1. Bloch, and S. Kuhr, Observation of Correlated Particle-Hole Pairs 
and String Order in Low-Dimensional Mott Insulators, Science 334 , 200 (2011); M. Endres, 
M. Cheneau, T. Fukuhara, C. Weitenberg, P. SchauB, C. Gross, L. Mazza, M. C. Bahuls, L. 
Pollet, 1. Bloch, and S. Kuhr, Single-site- and single-atom-resolved measurement of correlation 
functions, App. Phys. B 113 , 27 (2013). 

[67] Y. Ashida and M. Ueda, Diffraction-Unlimited Position Measurement of Ultracold Atoms in 
an Optical Lattice, Phys. Rev. Lett. 115 , 095301 (2015). 


26 



